Charge distribution in Ceo crystal doped by electric field 
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The calculations of the charge distribution in C6o-based 
FET structure are presented. A simple model is proposed 
to describe the distribution of the injected electrons or holes 
between two-dimensional layers. The calculations show that 
the relative charge distribution between the layers turns to 
be independent on the total amount of injected charges. The 
charge density is maximal on the surface layer and drops ex- 
ponentially with the depth increase. The relative amounts of 
injected charge involved in the top layer are 73 and 64 per 
cent in the case of electron and hole injection, respectively. 
Thus, the charge localization on the crystal surface turns to 
be markedly different from near complete one that was ob- 
tained earlier within tight- binding model for the charge con- 
centration providing superconductivity. 



I. INTRODUCTION 



is below 50 meV Thus, the fullerite crystal, which 
characterized by the relation U/t more than 20, have to 
be classified as a system with very strong electron cor- 
relation. That is why the results obtained with the use 
of the one-electron tight-binding model should be taken 
with care. 

The above reasons motivate the present study of 
charge-profile problem in fullerite. We consider the prob- 
lem from an opposite point of view that seems to be in 
better correspondence with the electronic characteristics 
of this material. We apply the approach earlier used by us 
for calculations of charge transfer in the fullerides A3C60 
crystals It is hopping smallness in comparison with 
electrostatic energies that is used to develop a simple 
scheme allowing us to calculate the charge distribution 
in the Cgo-based FET structure. 



Recently the superconductivity of organic substances 
due to injection of charges by electric field in the field 
effect transistor (FET) was found (]]-§]. This new phys- 
ical phenomenon is considered to be unusual in several 
aspects. First, unusual is the nature of the supercon- 
ducting objects themselves: they are organic crystals, 
anthracene, phenanthrcnc and fullerene Cgo- Second, un- 
usual is the superconductivity of the charges injected by 
electric field: up-to-date the superconductivity of organic 
crystals was related to the charge transfer, particularly, 
due to chemical doping. Third, in the case of fullerene 
Ceo the superconductivity is observed at extremely high 
temperature (52 K for pure Cgo @ and at even 117 K 
when doping by bromoform (^]). Fourth, injecting holes 
into the fullerene leads to significantly higher critical tem- 
perature than at injecting electrons, although at chemical 
doping the Cgo superconductivity is observed with nega- 
tive charge carriers only. 

These findings give rise to the problem of physical 
mechanisms of the phenomena, in the first turn, includ- 
ing the question on the distribution of charge carriers 
in organic crystal. In the study [|l| the amount of in- 
duced charges per Ceo molecule is estimated with the 
supposition that the additional charges occupy the sur- 
face monolayer only. 

An almost complete surface localization of induced 
charges near superconducting concentration is predicted 
in 0]. This study makes use of tight-binding model tak- 
ing into account electron repulsion by mean field. How- 
ever, in the case of fullerite, the Hubbard repulsion Uq 
on a Ceo molecule is estimated to be 1 eV or more || 
whereas the magnitude of the interfullerene hopping t 



II. MODEL AND METHOD 

As in Ref. [0, we consider the fee fullerite lattice with 
a (001) plane parallel to the gate. We are interesting 
in the distribution of injected charges between the (001) 
parallel layers of Ceo molecules that forms square lattices 
with the side length &=10 A. 

Let p n be total number of extra electrons per molecule 
in layer n, so their charge is —ep n (the case of injected 
holes is represented by negative p n ). We consider the 
system with the total numbers of layers N+l, the surface 
layer is labeled by zero. 

The basic assumption of our model is neglecting the 
overlap between Cgo molecules. Accordingly, the total 
energy of the crystal is the sum of the energies of the sep- 
arate layers of molecules that interact electrostatically. 
The energy E n per molecule in a layer n with the poten- 
tial U n is 



N 



E n = ^[E (pn) - eU n p n ] 



(1) 



where Eo(p n ) is energy of free molecule with p n extra 
electrons. 

This expression is accurate provided p n is integer, so 
Eo(pn) can be treated as energies of Cgo ions with corre- 
sponding charges. For noninteger p„, it is naturally as- 
sumed that Eo(p n ) can be interpolated by the quadratic 
fit, 
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E (p n ) = E (0) + E lPn + -E 2 p: 



(2) 
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Then the minimum condition of total energy with respect 
to p n leads to equation 



Ei + E 2 p n - eU n 



(3) 



where 9 is a Lagrange multiplier taking into account the 
additional restriction of fixing total electron number 



, Pn — Ptot- 



(4) 



From the other side, the layer potentials are, in turn, 
related to the charges by the Poisson equations. Because 
of negligible hopping, as it is adopted in our model, the 
Cgo molecules interact only electrostatically and can be 
represented by point charges due to their near spherical 
form. According to Ref . Q , the potential U n of layer n 
is determined by the equation 



eU n 



N 

£ ^2 min(m,n)p 

rn—0 



VPn, 



(5) 



where the coefficients £ and r\ are expressed through Cgo 
dielectric constant e and the distance b between neigh- 
boring molecules 
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(6) 



Excluding U n from (3) and (5) one obtains the system of 
linear inhomogeneous equations that determine, together 
with the condition (4), the charge densities p n , 



Xp„ 



m=0 
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U ^ Pn 

m— n+1 



a 



0,1. ..,7V, (7) 



where parameter A is defined by the equality 

x = (E 2 - n)/t 



(8) 



and C = (6 — E\)/£ is a new constant. As it can be seen 
from Eq.(7) at n = 0, this constant is in close relationship 
to the charge density in surface layer, C — Xpo- 

To treat the system (7), let us write it in the matrix 
form 



(A + AI)p = CV, 



(9) 



where I is unit matrix, V is the column vector all the 
components of which are units, and A is (N+1) x (N+1) 
matrix of form 



A = 
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(10) 



To find p n , let us solve the eigenvalue problem of the 
matrix A, 



X+AX = e, X+X = I, 



(11) 



where X is square matrix the column of which are eigen- 
vectors of A and e is diagonal matrix of corresponding 
eigevalues. Then one can obtain from (9) and (11) 



p = CX(e + AI) -1 X + V. 



(12) 



Eq.(12) shows that the shape of charge profile defined 
by the relative charge distribution, is independent on p tot . 
Obtaining p n at any C, they are normalized to ptot-, thus 
satisfying the additional condition (4). 



III. NUMERICAL RESULTS AND DISCUSSION 

The quantity 1 — p n / Ptot determines the relative part of 
injected charge below the surface layer of the crystal. Its 
dependence on the parameter E 2 is depicted in Fig.^, as 
the result of calculations with N=50. We consider only 
range E 2 > i] where the solutions have physical meaning 
and all the p n has the same sign as the constant C in (12). 
The Fig.0 shows that the charge density at (Ez—rj) — > +0 
is completely localized on the surface molecules. As E 2 
rises, the charges markedly expands into deeper layers. 
To obtain numerical estimations in real Cgo crystal one 
should evaluate only E 2 in the approximation (2). Con- 
sidering the latter as a fit to the ground-state energies 
of Cgo molecule and its ions one obtains the following 
estimations of E 2 for injected charges of both signs 



Eo 



Eo 



E(0) + E(2) - 2E(1) 

A(C 60 ) - A(Cq ) = 2.7eV, (for electrons) 

E(0) + E(-2) - 2E(-T) 

-I(C 60 ) + I(C+) = 3.8eVs (for holes) 



where A are electron affinities and I are ionizationpo- 
tentials with known experimental values (in eV) 



A(C 60 ) = 2.7, A(C, 
I(C 60 ) = 7.6, I(C, 
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7.6 = 11.4. 



Using the estimation of parameter £ in (6) £ = 2.9 eV [H, 
we obtain r\ = 3.9^^2/4^ =1.27 eV. With above estima- 
tions of E 2 the degree of charge derealization 1 — p n j p to t 
is evaluated as 0.266 for electrons and 0.359 for holes, as it 
is illustrated in Fig.[j]. In other words, about one quarter 
of injected electron charge and one third of hole charges 
turns out to be beneath the surface layer of Cgo crystal. 

The calculated charge distributions between the layers 
for electrons and holes in the TV = 50 crystal is shown 
in Figj2[ It is clear that the injected charge densities fall 
with the distance from the surface strictly exponentialy 
in both cases, with the electron density dropping more 
rapidly than hole one. 

To conclude, we summarize the main differences of our 
results from obtained in Ref. (i) The charge profile 
in m substantially depends on the charge concentration 
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whereas our relative distribution is independent of it. (ii) 
According to 98 % of the total charge is confined to 
the surface layer for doping p to t higher than 0.3, in op- 
posite to our calculations showing this degree less than 
75 % at any doping. The both results are in qualitative 
agreement with respect to domination of injected charges 
in the top layer. At the same time, our calculations show 
the possible nonnegligible occupation of undersurface lay- 
ers. The question on their role in the conductivity prop- 
erties of fullerite in FET structures needs in a separate 
consideration in future. However, just now our results 
suggest that the existing estimations of charge concen- 
tration providing the highest T c in superconductive FET 
structures may be subject to refinement. 



[1] J.H. Schon, Ch. Kloc, B. Batlogg, Nature 406, 704 (2000). 
[2] J.H. Schon, Ch. Kloc, B. Batlogg, Nature 408, 549 (2000). 
[3] J.H. Schon, Ch. Kloc and B. Batlogg, Science 293, 2432 
(2001). 

[4] S. Wehrli, D. Poilblanc, T.M. Rice. .Eur. Phys. J.B 23, 
345 (2001). 

[5] O. Gunnarsson, Rev. Mod. Phys. 69, 575 (1997). 

[6] V.A. Kupricvich, O.V. Shramko, O.L. Kapitanchuk, 

Mol.Cryst.Liq.Cryst. 361, 31 (2001). 
[7] H. Steger et al., Chem. Phys. Lett. 194, 452 (1992). 
[8] J. Hrusak and H. Schwarz, Chem. Phys. Lett. 205, 187 

(1993). 

[9] R.L. Hettich, R.N. Compton, R.H. Ritchie, Phys. Rev. 
Lett. 67, 1242 (1991). 



FIG. 1. Relative occupation of undersurface layers, 
1 — Pn/ ptot, vs curvature E2 of energy-charge fit (2). Dashed 
lines indicate the results of calculations for injected electrons 
and holes in C6o-based FET structure. 



FIG. 2. Charge density distributions (in logarithmic scale) 
calculated for injected electrons and holes in C6o-based FET 
structure. 
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